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ABSTRACT 

NGC 300 X-1 is a Wolf Rayet -I- black hole binary that exhibits periodic decreases in X-ray flux. We 
present two new observations of NGC 300 X-1 from the Chandra X-ray Observatory (totaling ~130 
ks) along with ACS imaging data from the Hubble Space Telescope. We observe significant short-term 
variability in the X-ray emission that is inconsistent with an occultation by the donor star, but is 
consistent with structure in the outer accretion disk or the wind of the donor star. We simultaneously 
fit a partially-covered disk blackbody and Comptonized corona model to the eclipse egress and non¬ 
eclipsing portions of the X-ray spectrum. We hnd that the only model parameters that varied between 
the eclipse egress and non-eclipsing portions of the spectra were the partial covering fraction (^86% 
during eclipse egress and ~44% during non-eclipse) and absorbing column ('^12.3x10^^ cm“^ during 
eclipse egress, compared to ^1.4x10^^ cm“^ during non-eclipse). The X-ray spectra are consistent 
with the movement of the X-ray source through the dense stellar winds of the companion star. From 
our new HST imaging, we find the WR star within the X-ray error circle, along with additional 
optical sources including an AGB star and an early-type main sequence star. Finally, we use our 
egress measurement to rephase previous radial velocity measurements reported in the literature, and 
find evidence that the velocities are strongly affected by the ionization of the wind by the compact 
object. Thus, we argue the inferred mass of the black hole may not be reliable. 

Subject headings: X-rays: binaries, individual (NGG 300 X-1) — accretion, accretion disks — stars: 
black holes 


1. INTRODUCTION 

NGC 300 X-1 (hereafter, X-1) is a high-mass X-ray 
binary (HMXB), l ocated 2.0 Mpc away in the spiral 
galaxy NGC 300 (|Dalcanton et al.l I2009D . The X-ray 
and optical emission has been studied by multiple au¬ 
thors, and all observations are consistent with the sys¬ 
tem cons isting of a black hole (BH) -b Wolf-Rayet (WRj 
binary (iCarpano et al.l l2007al: iCrowther et al.l 2?)?)^ 


Barnard et alJl2008l : lCrowther et al.1120101: [Binder et alJ 
20111) . The BH pri mary has been estima ted to have 


a mass of 20±4 Mq (|Crowther et al.ll2010D . making it 
the second largest stellar mass BH known (after another 
WR+ BH binary IC 10 X-1, iSilverman fc Filipp enkol 
I2008D . Heavy stellar mass BHs are expected to form 
in low metallicity environments, and both NGC 300 and 
1C 10 have measured metallicities that are ^ 15-30% solar 
(|Crowther et al.l[2003l : lUrbaneia et al.ll^OOSD . 

While numerous other BH-HMXBs with massive 
main sequence companions h ave been observed (see 
iMcClintock fc RemillardI 120001 for a review), BH -|- 
WR systems are exceedingly rare. Only one other X- 
ray source, IC 10 X-1, is a confirmed BH -|- WR binary 
(ISilverman fc FiliDpenkol 120081: iClark fc Crowtheil 12004 
iBauer fc Brandtll2004D . w hile Cyg X-3 in the M ilky Way 
is a BH -I- WR candidate (jLommen et al.l[2?)f)^ . Popula¬ 
tion synthesis models suggest that ~1 BH -|- WR binary 
may exi st as a bright X-ray source in a Milky Way-sized 
galaxy ( Ivan Ke rkwiik et al.l 119961 : iLommen et al.l 120051 : 
iZdziarski et al.ll2013 ). 

The orbital period of X-1 is short, ~33 hr 
(jCaroano et al.1 l2007bl : ICrowther et al.l 120101) . Periodic 


^ University of Washington, Department of Astronomy, Box 
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decreases in the X -ray flux, observed w ith both Swift 
and XMM-Newton (jCaroano et al.l[2007ar) . have been at¬ 
tributed to a grazing eclipse by the donor star. This 
has led to a rough constraint on the inclination angle of 
the system through geometrical argume nts (not formal 
dyna mical modeling) of z = 60 — 75° (jCrowther et al.l 
120101) . While similar periodic decreases in flux have been 
observed in other s ources (e.g., 1C 10 X-1 and other 
Galactic sources; see iBarnard et al.l[2?)T^ and references 
therein), it is not yet known whether the X-1 light curve 
is revealing an eclipse (a complete or partial occultation 
by the donor star) or a dip (where changes in the X-ray 
flux are due to structure in the donor star winds or an 
extended corona). Light curves of eclipsing sources are 
energy-independent, while dipping sources show periodic 
changes in X-ray hardness that correspond to changes in 
the flux level. 

We have obtained two new, deep Chandra observations 
of NGG 300, both of which contained X-1 in the field of 
view. Additionally, we have obtained optical imaging of 
X-1 using the Hubble Space Telescope {HST), allowing 
us to study both the X-ray and optical sources in detail. 
In this paper we present our observations and data re¬ 
duction techniques (jJH) and an analysis of the observed 
X-ray variability (jjS]) and spectra (jJH). We then discuss 
the implications for mass accretion onto the BH and the 
structure of the WR winds (j}5]), and conclude with a 
summary of our Hndings (jj6|). 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. X-ray Observations from Chandra 

We have obtained two new, deep observations of the 
spiral galaxy NGG 300 using the Chandra Advanced 
CCD Imaging Spectrometer (ACIS-I) taken on 2014 May 
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TABLE 1 

Chandra Observation Log & Alignment to USNO-Bl.O 


Obs. ID 

(1) 

Date 

(2) 

Effective Exposure 
Time (ks) 

(3) 

Read 

Mode 

(4) 

Data 

Mode 

(5) 

# Sources Used 
for Alignment 
(6) 

rms Residuals 
(arcsec) 

(7) 

Percentage 

Improvement 

(8) 

12238 

2010 Sept. 24 

63.0 

timed 

vfaint 

7 

07267 

55.44% 

16028 

2014 May 16-17 

63.9 

timed 

faint 

5 

07429 

7.01% 

16029 

2014 Nov. 17-18 

61.3 

timed 

faint 

9 

07423 

20.98% 


TABLE 2 

Summary of Chandra X-1 Observations 


Obs. ID 

(1) 

R.A. 

(J2000) 

(2) 

Deck 

(J2000) 

(3) 

Significance 

{^) 

(4) 

net counts 
(0.35-8 keV) 
(5) 

bkg. counts 
(0.35-8 keV) 
(6) 

Off-Axis 
Angle (°) 

12238 

00:55:10.00 

-37:42:12.2 

177 

2171±48 

125±1 

4.86 

16028 

00:55:10.02 

-37:42:12.3 

273 

2501±51 

66±1 

4.85 

16029 

00:55:9.99 

-37:42:12.1 

225 

1692±42 

41±1 

4.84 


16-17 and 2014 November 17-18. X-1 was contained 
within the field of view of both observations. We sup¬ 
plement our new observations with a publicly available 
Chandra AClS-1 observation of NGC 300 from 2010 
(|Binder et al.l 1201111 . All X-ray observations were re¬ 
duced using the Chandra Interactive Analysis of Ob¬ 
servations (CIAO) software package version 4.6.1 and 
CALDB version 4.6.3 using standard reduction proce¬ 
dures. The data from all three observations were re¬ 
processed from evtl using chandrajrepro. Background 
light curves were extracted and no strong background 
flares were found. We used the task lc_clean and 5(t 
clipping to create good time intervals (GTIs) for each 
observation; all event data were filtered on the GTIs. 
Table [I] provides an observation log of all three observa¬ 
tions, including: the observation ID number (hereafter 
referred to as the “Obs. ID”), the date of the observa¬ 
tion, the effective exposure time in kiloseconds, and both 
the read mode and data mode. 

The evt2 data were aligned to the USNO-Bl.O cat¬ 
alog using the CIAO tasks wcsjnatch and wcs_update. 
Exposure maps were constructed using the CIAO tool 
fluximage, which automatically produces exposure- 
corrected images using a user-specified instrument map. 
We assume spectral weights appropriate for both XRBs 
and AGN: a power-law spectrum (with T = 1.7) ab¬ 
sorbed by the ave rage foreground colu mn density (A^h = 
4.09x10^° cm~^. iKalberla et al.l[200!Tti . We additionally 
use the task reproject_image to merge all three expo¬ 
sures into a single image. Point sources in each exposure 
and the merged image were detected using wavdetect, 
and sources with a signal-to-noise ratio > 3 were consid¬ 
ered significant. 

X-1 was detected at high significance in all three ob¬ 
servations, at positions in excellent agreement with the 
position reported in SIMBAEQ. Table [2] summarizes the 
observed properties of X-1 in each individual exposure, 
including the wavdetect R.A. and Deck, the off-axis an¬ 
gle from the nominal aim-point in each observation, de¬ 
tection significance, and net and background counts (in 
the 0.35-8 keV band). We use PIMM^ v4.2 to estimate 
the level of pile-up suffered during each observation and 

^ See http://sinibad.u-strasbg.fr/ 

® See http: //cxc. harvard. edu/toolkit/piimns. j sp 


find the fraction to be low, < 5%. 


Fig. 1.— Chandra 0.35-8 keV image of X-1 in ObsID 16029. The 
red cross indicates the position of the X-ray source, while the white 
circle shows our source extraction region and the yellow annulus 
shows the region used for background extractions. 
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Fig. 2.— The radial surface brightness profile of the merged X- 
1 image. The dashed gray line shows the background level, and 
the vertical dotted line shows the inner radius of the background 
extraction region. 


Figure [T] shows a 0.35-8 keV image of X-1 from Ob- 
sID16029 (the observation from ObsID 16028 looked sim¬ 
ilar). We find no evidence for diffuse or extended emis¬ 
sion in the vicinity of X-1. The source region shown 
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in white (used to extract light curves and spectra) was 
chosen to be a circular region, centered on the merged 
source coordinates, with a radius that enclosed ^90% 
of the Chandra PSF at that position. To determine the 
background region (the circular annulus shown in yel¬ 
low), we first extracted a radial surface brightness profile 
from the merged Chandra image of X-1. The inner radius 
was set at the distance from the source where the radial 
surface brightness profile reached background levels, as 
shown in Figure[2l The outer radius was then determined 
such that the source and background regions contained 
approximately the same number of pixels. Source and 
background regions were examined by eye to ensure that 
no other obvious sources were observed in either region. 


2.2. Optical Observations from HST 

In addition to our two new Chandra exposures, we have 
new Hubble Space Telescope {HST) imaging to study the 
resolved stellar populations in the immediate vicinity of 
X-1. To directly compare our X-ray observations with 
the optical ones, we first needed to place both sets of 
observations on the same coordinate system. 

Directly aligning both the Chandra and HST images 
was not reliable due to the difficulty matching X-ray 
sources with optical counterparts in the HST field. We 
therefore performed relative astrometry by aligning both 
the X-ray and optical images to the same large-held 
optical reference image and coordinate system. To do 
this, we retrieved a publi cly-available S-band image from 
NEE0 (iKim et al.l l2004f) . The reference S-band image 
was aligned using the USNO-Bl.O catalog using the posi¬ 
tions of six stars. The IRAF task ccmap was used to com¬ 
pute the plate solution and to update the image header 
with corrected WCS information; the root-mean square 
(rms) residuals of the ht were 0"158 in right ascension 
and 0"159 in declination. We were next able to align our 
HST helds to the ground-based R-band image by identi¬ 
fying hve matched sources. The resulting rms residuals 
from the B-hand/HST alignment were 0"125 in right as¬ 
cension and 0"0625 in declination. Adding in quadrature 
the rms residuals for both helds produces a hnal optical 
astrometry rms of 0"201 in right ascension and 0"171 in 
declination. 

We searched for optical HST counterparts in the vicin¬ 
ity of the X-1. The error circle had a radius of 0'.'69; 
fourteen optical sources were found within the X-ray er¬ 
ror circle, including a likely WR star. The likely WR 
star is detected at an RA (J2000) of 00:55:9.99 and Deck 
(J2000) -37:42:12.65. This position is ~075 off from the 
previously reported position of the WR from ground- 
based observations ( Schi ld et al.l 120031: ICaroano et al.l 
120061 : iCrowther et al.1 1200711 . However, the seeing asso¬ 
ciated with the ground-based observations was on the 
order of ~1", and the region is crowded with other rel¬ 
atively bright stars, making our observed position con¬ 
sistent with earlier observations. In Figure [3] we show 
the position of the X-ray source and the WR star, with 
their associated error circles, for both the HST image 
and the ground-based S-band image. To create an RGB- 
rendered image from our HST observations, the FilAW 
image is set to the red channel, the FdOQW filter is used 


http://ned.ipac.caltech.edu/ 


for the green channel, and the blue channel is estimated 
by taking 2 x F606f^ - F814Vk. 



Fig. 3. — Left: Smoothed, 5^^ x 5^^ ground-base S-band image 
of NGC 300 X-1. The position of the X-ray source and the cor¬ 
responding X-ray error circle are shown in blue, and position of 
the WR star is shown in red. Right: The optical, RGB-rendered 
HST image of NGC 300 X-1. The FSIAW exposure is shown in 
red, F606W is shown in green, and the blue band is estimated as 
2xG — R. The box is 5^' on a side. The position of the X-ray source 
and the corresponding X-ray error circle are shown in white; the 
WR candidate position and optical error circle are shown in red. 
The red arrow indicates the AGB star, and the cyan arrow indicates 
the likely high-mass main sequence star. 


To measure resolved stellar photometry on the HST 
imaging data, we used the same techniques and software 
as those applied to the ACS im aging of the Panchrom atic 
Hubble Andromeda Treasury (|Williams et al.l [201411 . In 
short, we used the point spread function photometry 
package DOLP HOT - an updated version of HSTPHOT 
(iDolnhinl 1200011 . In DOLPHOT, all of the individual 
CCD exposures are aligned and stacked in memory to 
search for any significant peaks, and then each signifi¬ 
cant peak above the background level is fitted with the 
appropriate point spread function. The measurements 
are corrected for charge transfer efficiency and calibrated 
to infinite aperture. The measurements are then com¬ 
bined into a final measurement of the photometry of the 
star. The final measurements include a combined value 
for all data in each band for the count rate, rate error, 
VEGA magnitude and error, background, x of tho PSF 
fit, sharpness, roundness, crowding, and signal-to-noise. 
For our photometry reported here, we use the combined 
VEGA magnitude in each observed band. 

We measure the F606W and FSIAW magnitudes of 
the WR candidate to be meoe = 22.412 ± 0.005 and 
771814 = 22.327 ± 0.007, consistent with ground-based 
esti mates of the Vyand a pparent magnit ude of 22.53 
from iCarpano et al.l (|2006f) and 22.44 from iSchild et al.l 
( 200311. Using a diy ance of 2 Mpc to NGC 300 
(|Dalcanton et al.l[2009ll . these apparent magnitudes cor¬ 
respond to absolute magnitudes in FQOCW and F9>1AW 
of -4.13 and -4.20, respectively, assuming Galactic fore- 
grou nd extinction values of Av= 0.035 and A/=0.019 
fromiSch laflv &: Einkbeineil (|2011ll , a recalibration of the 
(iSchlegel et al.l 1199811 dust maps. Optical spectroscopy 
of this star from the VLT/FORS 2 shows a WR stm wiH i 
an early-type WN 5 spectrum (jCrowther et al.l 1200711 , 
which typically have a mean intrinsic absolute V mag¬ 
nitude of_;T2_J^^£ls|^Tiffl£311003) . Erom the spec¬ 
trum, ICrowther et al.l (|2010ll estimate the surface tem¬ 
perature of the star to be T* =65000 K, with a luminosity 
log(L/L0)=5.92 and a wind speed 7;oo=1300 km s“^. 

In addition to the WR candidate, thirteen other opti- 
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cal sources are detected within the X-ray error circle. A 
color magnitude diagram (CMD) of all the sources de¬ 
tected within the HST field of view is shown in Figured 
with the sources falling within the X-ray error circle su¬ 
perimposed in purple. The WR candidate is shown in 
blue, and is clearly brighter and bluer than the other 
sources found within the error circle. Another early-type 
star found in the X-ray error circle is ~2.5 mag fainter, 
and likely a main sequence star (late B-type based on an 
FQQQW-FSIAW color of ~0). Another very bright star 
(likely an AGB) is also consistent with the X-ray posi¬ 
tion. It is therefore possible that the WR candidate is 
not the true optical counterpart to the X-ray source. In 
(J5]we discuss the possibility that the optical companion 
to the X-ray source is a massive main sequence star or 
an AGB star. 

X-ray-to-optical flux ratios, ^og{fx/fv), can be a 
useful tool for discriminating between different classes 
of astrophysical objects; for example, in the SMG 
the X-ray-to-optical flux ratio combined with an op¬ 
tical c olor separates NS-HM XBs from LMXBs and 
AGN (|McGowan et al.l I2008D . We therefore compute 
log(/x//y) for each optical counterpart candidate to in¬ 
vestigate whether any of the optical counterpart candi¬ 
dates are consistent with a NS-HMXB origin. We use the 
F606W band and 2-10 keV band fluxes as fv and fx, 
respectively. Table |3] lists the distance from the X-ray 
source position, the optical magnitudes, and ^og{fx/fv) 
values, assuming that optical source were the true coun¬ 
terpart to the X-ray source. 

The NS-HMXBs in the Small Magellanic Gloud all 
exhib it X-ray-to-optical flux ratios <1 (|McGowan et al.l 
1200811 . The X-ray-to-optical flux ratios computed here 
provide further evidence that, even if the WR star is not 
the true companion to the X-ray source, X-1 is inconsis¬ 
tent with a NS-HMXB origin. However, with the large 
values of ^og{fx/fv) and red colors, many of the possi¬ 
ble counterparts would be consistent with an LMXB ori¬ 
gin. Thus, we cannot rule out the possibility of a LMXB 
origin from the optical photometry alone. We further 
discuss the nature of the optical counterpart in H5.2I 


TABLE 3 

Properties of Optical Counterpart Candidates 


Distance from 
X-ray source 
(1) 

UIF606IV 

(2) 

(3) 

log(/x//F 6061 v) 

(4) 

o'.'451 

22.412±0.005 

22.327±0.007 

1.38 

O'.'349 

23.726±0.019 

22.109±0.014 

1.91 

O'.'108 

24.948±0.024 

23.673±0.015 

2.40 

O'.'075 

25.049±0.026 

24.100±0.024 

2.44 

0'.'127 

24.808±0.022 

24.757±0.038 

2.34 

0'.'519 

25.558±0.032 

24.669±0.029 

2.64 

0'.'053 

26.183±0.049 

25.707±0.062 

2.89 

0'.'260 

26.515±0.066 

25.449±0.051 

3.03 

0'.'164 

26.580±0.089 

25.399±0.065 

3.05 

0'.'051 

26.248±0.071 

25.647±0.087 

2.92 

0'.'049 

26.739±0.106 

26.051±0.118 

3.12 

(X'249 

26.678±0.093 

26.068±0.109 

3.09 

0'.'136 

27.178±0.156 

26.202±0.135 

3.29 

O'.'224 

27.681±0.187 

26.967±0.191 

3.49 


Note. — The first row is the likely WR star. 
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F606W-F814W 

Fig. 4. — A color-magnitude diagram showing all the stars within 
the HST field. Purple circles indicate the optical sources found 
within the Chandra error circle. The point representing the WR 
star is shown in blue, and is significantly brighter and bluer than 
the other stars within the error circle. The red arrow indicates the 
position of the AGB star, and the cyan arrow indicates the position 
of the likely high-mass main sequence star. 

3. X-RAY VARIABILITY 

Variability is a common feature of many HMXBs, and 
the timescales over which a single source changes its X- 
ray luminosity can range from minutes to decades. X-1 
has been previously observed to be variable, with dips in 
the X-ray light cur ve attributed to a glanc ing eclipse with 
a period of ^^33 hr (|GarDano et ahll^OTall . Geometric ar¬ 
guments have been used to constrain the most likely in¬ 
clination of the system to be z = 60 — 75° (|GarDano et al.l 
I2007at IGrowther et al.l[2010ll . 

The 0.35-8 keV light curves from our Chandra obser¬ 
vations were extracted using the GIAO tool dmextract. 
The cumulative photon arrival time distributions are 
shown in Figure [51 The time axis is shown as a per¬ 
centage of the total exposure time to account for the 
slight differences between the observations. The light 
curves, binned to 3 ks, are shown in Figure El In ObsID 
16028, X-1 showed a secular increase in 0.35-8 keV count 
rate. A two-sided KS test against a constant count rate 
(equal to the mean count rate in the observation) yielded 
a probability of a constant count rate of 7.5x10“^%. In 
ObsID 16029, however, the observed light curve is con¬ 
sistent with a constant count rate and looks similar to 
the light curve observed in ObsID 12238. 

3.1. Periodicity 

The three light curves were folded over a variety of 
orbital periods allowed by the earlier Swift and XMM- 
Newton observations (a range of 32.4-33.2 hr, in steps of 
0.1 hr). The start of ObsID 16028 was set to phase zero 
and matches closely wi th the secular increase observed in 
(iGaroano et al.l[2007aL their Figure 4). Figure [71 shows 
the light curve folded on a period of 33.0 hr; folding on 
different periods produced similar results. The first ^20 
ks of ObsID 12238 may represent eclipse egress, although 
the majority of the ObsID 12238 and 16029 exposures 
were taken during phases approaching conjunction. 

To determine whether X-1 exhibited evidence of pe¬ 
riodic X-ray pulsations on timescales shorter than the 
exposure time, we performed a periodogram analysis 
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Fig. 5.— The cumulative 0.35-8 keV photon arrival time distri¬ 
butions for both Chandra observations. The dashed line shows the 
expected distribution for a constant count rate. The colors indi¬ 
cate the observation ID: 12238 is shown in blue, 16028 is shown in 
purple, and 16029 is shown in red. 

(|Home fc Baliunaslll986D of each exposure. Monte Carlo 
simulations were used to calculate the 68%, 90%, and 
99% confidence levels assuming the light cu rve was dom¬ 
inated purely by Poisson noise. While [Binder et al.l 
(|2011ll reported one peak at 0.785 hours that exceeded 
the 90% confidence level (but was determined unlikely 
to be a genuine periodic signal), we find no obvious 
peaks in the X-1 periodograms above the 90% confidence 
level. The power density spectrum does not show signif¬ 
icant differences from previous observations, and is well- 
described by a power law with a spectral index of ~1. 
This spectral index is characteristic of accr eting BHs in 
eithe r the thermal or stee p power law state (Ivan der Klisl 
11994 [19^ iBellonil 120101 1. We additionally searched for 
periodicity in the hardness ratio light curves, but did not 
find any significant peaks. 

3.2. Spectral Variability 

The hardness ratio light curves for each observation 
are additionally shown in Figure |6l If the dip in the 
X-1 light curve observed in ObsID 16028 was due to 
the donor star partially eclipsing the BH, there would 
be no change in the observed X-ray spectrum; one 
would simply observe the same spectral shape at a lower 
flux level. Any energy-dependent evolution in concert 
with the light curve would suggest the X-ray source 
is not being eclipsed by the companion star, but by 
dense stellar winds or structure in the outer accretion 
disk. Similar beh avior has been observed in Galac¬ 
tic eclipsing XRBs |Church fc Balucihska-Churchl [20011 : 
iBarnard et al.l 1200111 anT more recently, IG 10 X-1 
(jBarnard et al.lT2ni4 l. To search for energy-dependent 
evolution, we extracted light curves in the soft 0.35-2 
keV and hard 2-8 keV energy bands. We then compared 
the ratio of hard/soft band counts to the total 0.35-8 keV 
count rate. We found no evidence of a change in hardness 
in ObsIDs 12238 and 16029, but ObsID 16028 exhibited 
an elevated hardness during the first ^^20 ks of the ob¬ 
servation, corresponding to the portion of the exposure 
when the sharpest increase in count rate was observed. 
The null hypothesis that the ObsID 16028 light curve is 
consistent with a constant hardness ratio is rejected at 


the 3.6(7 level. A more detailed investigation of the spec¬ 
tral evolution during ObsID 16028 is discussed in the 
next section. 

The hardness ratio peaks at the deepest point dur¬ 
ing the presumed eclipse egress and steadily decreases, 
reaching a relatively constant minimum level at superior 
conjunction with the donor star. Although our Chandra 
observations do not contain the eclipse ingress, we expect 
that additional observations during this orbital phase 
would reveal air increase in hardness ratio. At phases 
consistent with superior conjunction, there appears to 
be some low level of aperiodic variability, likely due to 
the BH moving through regions of the donor Wolf-Rayet 
stellar wind of changing optical depth; this scenario is 
discussed in more detail in |j5l 

4. SPECTRAL MODELING 

We extracted spectra and response files for both new 
X-I Chandra observations in the 0.35-8 keV band using 
the tool specextract and grouped to con tain at least 20 
counts per energy bin. We use XSPEC (|Arnaudl 1199611 
v.l2.6.0q to perform all spectral fitting; the goodness 
of each ht is evaluated using statistics and standard 
weighting. All errors correspond to the 90% confidence 
level, assuming symmetrical errors (e.g., = 2.71, 

4.61, and 6. 25 for one, two , or t hree free parameters, 
respectively; iLamoton et al.l 1197611 . We assume a dis- 
tance to NGG 300 of 2.00±0.04 Mpc (jPalcanton et al.l 
1200911 . and all models include neutral absorption due 
to the Galactic col umn of Ah = 4.09x10^° cm“^ 
(|Kalberla et al.l[2005ll . 

Previous studies with X MM-Newton (|CarDano et al.l 
l2007bl : IBarnard et al.l[2?i0^ have shown evidence that X- 
1 may undergo state transitions. In some observations, 
the X-ray spectrum is described by a power law with a 
photon index of ^4 contributing ^60% of the total X-ray 
flux; the remainder of the flux originates in a multicolor 
disk blackbody (with fcPin ^2 keV). This is consistent 
with a “thermal” BH state (jMcGlintock fc RemillardI 
l2006fl . However, in other observations the X-ray emis¬ 
sion resembles a steep power law (SPL) state, with a 
photon index near 2.5 and an emission line near 1 keV. 
A longer exposure taken with Chandra revealed a power 
law spectrum (P ~2), comprising ~74% of the 0.35-8 
keV flux, and a significantl y cooler disk blackb ody with 
a temperature of ^0.2 keV (|Binder et al.ll20fll l. 

Since the spectrum of X-1 has previously required both 
power law and thermal components, we focus on three 
models of the 0.35-8 keV spectrum: a single power law 
model, a combined power law with a disk blackbody, and 
a power law with some emission originating in a thermal 
plasma (apec). We choose to use the apec component 
because the model, when combined with a power law, can 
produce excess emission near 1 keV without the need for 
a Gaussian component. 

The results of our spectral fitting are shown in Ta¬ 
ble m Although a simple power law model yielded a 
X^/dof=103/98 for the spectrum from ObsID 16028, sys¬ 
tematic residuals were present near keV. The addi¬ 
tion of an apec component improved the fit residuals, 
with r =2.29±0.12 and kT = 1.3lg3 X^/dof = 

82/96. An /-test was used to determine if the addition of 
the thermal component provided a significant improve- 
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Fig. 6 .— Top: the 0.35-8 keV light curves of X-1, binned to 3 ks. The jdoj value each panel shows the result from fitting to a constant 
count rate. Bottom: the ratio of counts observed in the hard 2-8 keV band to counts observed in the 0.35-2 keV band, also binned to 3 
ks. ObsID 12238 is shown in blue (left column), 16028 is shown in purple (middle column), and 16029 is shown in red (right column). The 
solid horizontal line indicates the median count rates or hardness during each observation. Obs ID 16028 shows a clearly secular increase 
in the count rate, with a significantly harder spectrum during the first -^^20 ks of the observation. 
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TABLE 4 

X-1 Spectral Models^ 


Parameter 

(1) 

po 

( 2 ) 

po+diskbb 

(3) 

po+apec 

(4) 


ObsID 16028 



r 

fcTi„ (keV) 
kT (keV) 

% PL 

X^/dof 

flux*’ (erg cm~^) 

ftest probability against po 

2.39±0.07 

100 

103/98 

(5.8±0.1)x10-13 

2.45±0.40 

0 9+d'7 

92 

101/96 

(5.7±0.6)x10-13 

47% 

2.29±0.12 

. O+0.4 
-^••^-0.3 

92 

82/96 

(5.5±0.1)x10-13 

0.08% 


ObsID 16029 



r 

Q QQ“t"D.uy 

Z. JO_Q Q8 

'4-0.73 

9 9 C-I-U .10 

16 

fcTin (keV) 


1 4+0-6 
-*-•^-0.3 


kT (keV) 



1.1±0.3 

% PL 

100 

79 

94 

X^/dof 

66/68 

64/66 

57/66 

fiux^ (erg cm~^) 

(5.3±0.1)x10-13 

(5.5/“;|)xl0-« 

(5.0±0.2)xl0-13 

ftest probability against po 


36% 

0.8% 


Note. — "^None of our spectral fits required absorption beyond the Galactic column. We therefore kept A^h fixed at 4.09x10^^ cm ^ 
for all models. ^All fluxes are unabsorbed and measured in the 0.35-8 keV energy band. 

of NGC 300. Roughly 30% of the flux of X-1 originated 
in the power law component. 

The spectrum from ObsID 16029, however, was well 
described by a simple power law with T =2.38lQg3 
(x^/dof = 66/68). We found no evidence that adding 
either a disk blackbody or a thermal component was 
required by the observed spectrum. The best-fit 
spectrum yielded an unabsorbed 0.35-8 keV flux of 
(5.3±0.1)X10“^^ erg s“^ cin“^, or (5.2±0.1)x 10^® erg 


ment over a simple power law fit. The addition of the 
thermal component produces an /-statistic of 0.08%, in¬ 
dicating the addition of the third component improved 
the fit compared to the simple power law model. This 
model is consistent with the power law plu s Gau ssian 
emission line model ut ilized in iBarnard et al.l (j2008f ) and 
iCarnano et al.l (|2007bf ). and yielded an unabsorbed 0.35- 
8 keV flux of (5.5 ± 0.1) X10“^® erg s“^ cm“^, equivalent 
to a luminosity of (2.6±0.1)xl0®® erg s“^ at the distance 
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Fig. 7. — The 0.35-8 keV light curve, folded on a period of 33.0 
hr, and binned to 3 ks for clarity. This period is consistent with 
previous measurements from Swift and XMM-Newton. 

at the distance of NGC 300. 

The best-fit parameters for our various spectral models 
are summarized in Table |4l Our favored spectral model 
is printed in boldface. Figure [8] shows the 0.35-8 keV 
spectrum of X-1 in each observation (ObsID 16028 in 
the top row; ObsID 16029 in the bottom row), with ei¬ 
ther a simple power law model (left column) or a power 
law plus thermal emission model (right column) super¬ 
imposed. The fit residuals are shown immediately below 
each spectrum and model. 

To test whether our current data are consistent with 
previous models of the X- 1 spectrum, we used the 
best-fit mod e ls rep orted in iBarnard et al.l (|2008[1 and 
[Binder et al.l (1201111 . The fit parameters were frozen to 
their best-fit values, so that the only parameters allowed 
to vary were the respective normalizations. The resulting 
values are reported in Tabled ObsID 16028 was not 
well described by any of the previously reported best-fit 
models, while the ObsID 16029 spectrum showed good 
agreement with several earlier models. This lends fur¬ 
ther support to the notion that ObsID 16029 was ob¬ 
served during a similar orbital phase as Obs ID 12238 
and observations 2-4 of the XMM-Newton observations. 
In fact, we obtained the worst value when we mod¬ 
elled the ObsID 16029 spectrum with the be st-fit model 
from observation 1 of IBarnard et al.l (I2008I1 . when X-1 
was observed during eclipse and egress. 

4.1. Temporal Variations in the Spectrum of ObsID 

16028 

Variations in the X-1 ligh t curve are consistent with 
an orbital period of ^33 hr (jGarpano et al.l[2007ari . Al¬ 
though not fully captured with our Chandra observa¬ 
tions, the X-ray eclipse data from Swift and XMM- 
Newton indicate an eclipse duration of ^6 hr. Both the 
orbital period and eclipse duration are similar to those 
of IC 10 X-1, another eclipsing WR -I- BH binary, which 
additionally exhibits spec tral changes during eclipsing 
and non-eclipsing periods (iStrohmaver fc Pashamll2013l : 
IBarnard et al.l[^14l : iLavcock et al.ll2015h . The evolution 
of the IC 10 X-1 spectrum with orbital phase provides 
strong evidenc e that the system ho sts a substantial ex¬ 
tended corona (|Barnard et al.l[201^ . 

The sharp increase in count rate observed during the 


first ^20 ks (^5.5 hours) of ObsID 16028 is consistent 
with previous observations of the X-1 eclipse egress. Fur¬ 
thermore, the observed energy evolution in the light 
curve indicates that we are not observing a simple oc- 
cultation of the X-ray source by the photosphere of the 
companion star, and that other structures in the binary 
(stellar winds, an extended corona, etc.) are likely con¬ 
tributing to the observed periodic decrease in flux. 

To investigate the cause of the changing hardness as a 
function of time, we extracted two spectra from the Ob¬ 
sID 16028 observation (during the eclipse egress, from 
0-20 ks, and from 40 ks until the end of the observation, 
when the X-ray source was roughly halfway through its 
orbit) to search for variations in the spectral shape with 
time. The egress spectrum was binned to contain at least 
10 counts per bin, as the first 20 ks contained only ^380 
net counts, while the final third of the observation con¬ 
tained a sufficient number of counts to allow the spec¬ 
trum to be binned to at least 20 counts per bin. 

We expect that the BH in X-1 hosts both an accretion 
disk, emitting as a disk blackbody, and a corona. We first 
attempted to determine whether the eclipse egress spec¬ 
trum was consistent with a compact or extended corona 
scenario. We used the XSPEC models diskbb (to model 
the X-1 accretion disk) and comptt (to model a Comp- 
tonized corona), with both components being subject to 
the same absorbing column. In the XMM-Newton ob¬ 
servations, the best-fit spectral models yielded high in¬ 
ner disk temperatures (^2 keV) and power law indices 
F > 3. In a compact corona scenario, such parame¬ 
ters are unphysical, as a compact coron a can only ac¬ 
cess photons from the inner acc retion disk ([Roberts et al.l 
l2005t iGoncalves k. Serial 1200611 . In our spectral model, 
the compact corona is represented by forcing the inner 
disk temperature to be tied to the seed photon temper¬ 
ature. In an extended corona scenario, the corona is 
large enough to access soft photons from the outer edge 
of the disk in addition t o the hot inner-disk photons 
(|Haardt fc Maraschilll993ll . therefore the disk tempera¬ 
ture and seed photon temperature are allowed to vary 
independently from one another. 

We first considered the eclipse egress portion of the 
spectrum. For the compact corona scenario, we find 
a best-fit inner disk temperature (and, therefore, seed 
photon temperature) of 0.16±0.11 keV, with 38% of the 
flux originating from the compact corona. The result¬ 
ing x^/dof was 24/25. When we allow the inner disk 
temperature and the seed photon temperature to vary 
independently (e.g., the extended corona scenario), the 
results do not change significantly: fcTin = 0.20lQ j!)g keV 
and the resulting y^/dof was 24/24. We note that we 
recover these fit parameters even when we initialize the 
fit with much hotter inner disk temperatures and lower 
seed photon temperatures. 

The low inner-disk temperatures derived in these mod¬ 
els imply a str ong retrograde black hole spi n, which 
is unlikely (see iMcClintock fc RemillardI 120061 and ref¬ 
erences therein). With the relatively low number of 
counts present in the egress spectrum (^380), we can 
not statistically differentiate between the compact and 
extended corona scenarios using the spectral model de¬ 
scribed above. As described in the previous section, how¬ 
ever, we do obtain thermal temperatures of '^1 keV when 
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Fig. 8 .— The 0.35-8 keV spectrum of X-1 during ObsID 16028 {top row) and ObsID 16029 {bottom row). The left column shows the 
spectrum with a simple power law model superimposed, while the right column shows a power law plus thermal emission model. 


TABLE 5 

Comparing the Current Data to Past Spectral Models: x^/dof 



Barnard et al. {'2008') 


Binder et al. ('2011') 

Obs. 1 

Obs. 2 

Obs. 3 

Obs. 4 

(1) 

(2) 

(3) 

(4) 

(5) 

ObsID 16028 139/90 

ObsID 16029 71/68 

103/91 

67/69 

114/91 

66/69 

119/91 

69/69 

112/90 

70/68 


using simpler spectral models and the full, time-averaged 
spectrum. We therefore expect that further observations 
of X-1 during the X-ray eclipse would show a prefer¬ 
ence for the extended corona scenario, as was observe d 
in IC 10 X-1 (|Barnard et al.l[2f)T^ iLavcock et al.ll2015H . 

We next compared the eclipse egress spectrum with 
the spectrum extracted from the final third of the Ob¬ 
sID 16028 observation. The X-ray eclipse may be due 
to the BH’s motion through the dense stellar winds of 
the WR star, with some partial occultation by the WR 
star. To test this scenario, both spectra were modelled 
as an accretion disk (diskbb, initialized with an inner 
disk temperature near 2 keV) and an extended corona 
(comptt). Both components were subject to the same 
absorbing column, and we applied an additional partial 
covering model to the corona. As when testing for dif¬ 
ferent corona types, assuming a compact corona (instead 
of an extended corona) did not significantly affect the 
quality of the resulting fits. 


Both the eclipse egress and non-eclipsing spectra were 
simultaneously fit with this model. All parameters (ex¬ 
cept the normalizations) for the disk, corona, and absorb¬ 
ing column were tied between the eclipse egress and non¬ 
eclipsing spectra. The covering fraction and column den¬ 
sity of obscuring material in our partial covering model 
component were allowed to vary. The best-fit model is 
shown in Figurejlland had y^/dof of 69/70. The eclipse 
egress spectrum is shown in black and the non-eclipsing 
portion of the spectrum is shown in red. The best-fit in¬ 
ner disk temperature kT-m = 1-1 ±0.9 keV and the seed 
photon temperature was 0.1 keV. The electron tempera¬ 
ture was ^40 keV and the optical depth was 0.17. Nei¬ 
ther spectra showed evidence of requiring an absorbing 
column beyond the Galactic column, and the unabsorbed 
fluxes of the disk and corona did not change: the total 
unabsorbed 0.35-8 keV flux was ~ 9.6 x 10“^^ erg s“^ 
cm“^, corresponding to a luminosity of 4.6x10^® erg s“^, 
with ~90% originating from the Comptonized corona. 
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Fig. 9.— The 0.35-8 keV eclipse egress spectrum (black) and 
the non-eclipsing spectrum (red). Both spectra were fit with a 
disk blackbody and a partially covered Comptonized component 
representing the extended corona. The only parameters allowed 
to vary between the two fits were the relative fiuxes of the disk 
and corona, the column density of the absorber, and the covering 
fraction. The corona was suffered more absorption and a higher 
covering fraction during the eclipse egress spectrum as compared 
to the non-eclipse spectrum. 

The partial covering component of our spectral model 
exhibited variability between the two spectra. During the 
eclipse egress, the partial covering fraction was ( 86 ^ 21 )% 
with an absorbing column of iVH,cov=( 12 . 3^2 g) ^ 10 ^^ 
cm“^. During the non-eclipsing portion of the spectrum, 
the covering fraction dropped nearly in half ( 44 ^ 15 )%, 
while the absorbing column dropped by almost an or¬ 
der of magnitude to (1.4l^ g) x 10^^ cm“^. These re¬ 
sults (summarized in Table [H) suggest that the observed 
change in flux is due to a cloud of absorbing material 
periodically obscuring most of the X-ray source during 
the binary’s orbital period. We discuss the X-1 corona 
further in i l5.ll 

5. DISCUSSION 
5.1. An Extended Corona? 

Our spectral modeling provides preliminary evidence 
that X-1 may host an extended, hot corona. The tem¬ 
perature of the inner edge of the accretion disk, column 
density of obscuring material, and inferred luminosities 
are all remarkably simila r to those observed for IC 10 
X-1 (|Barnard et al.l 120141) . While our Chandra observa¬ 
tions of X-1 are not deep enough to definitely rule out the 
compact corona scenario, the many similarities to IC 10 
X-1 (which can be observed in much greater detail due 
to its closer distance) suggest that the X-ray spectrum 
of X-1 is being produced in a similar fashion. 

Observations by Chandra and XMM-Newton of IC 10 
X-1 and other dipping XRBs often reveal asymmetric 
eclipses, with the eclipse ingress being steeper (spanning 
a shorter period of time) than the egress. While our 
Chandra observations do not span a full orbital cycle, we 
can estimate the duration of the eclipse ingress and egress 
from the X MM-Newton and Swift observations of X-1 
presented in iCarpano et al.l (l2007all : the deepest portion 
of the eclipse spans ^ 0.2 in orbital phase (^^ 6.6 hours), 
while the ingress and egress occur over a phase period 


of ~0.25 (^8 hours) and ~0.15 (^5 hours), respectively. 
Unlike other dipping sources, the NCC 300 X-1 ingress 
is longer than the egress. 

Usin g the approach of iChurch fc Bahicihska-Churchl 
(l 200 l . we estimate the size of the extended corona from 
the eclipse ingres s using the estimatio n 27rrDA< = 2rcP 
(e.g., equation l. lBarnard et al.ir2014D . In this equation, 
P is the orbital period. At is the eclipse ingress time, 
CD is the radius of the accretion disk, and rr, is t he ra - 
dius of the corona. IChurch fc Balucihska-Churchl (j2004f ) 
assume the disk radius cd to b e 80% of the Roche lobe 
equivalent radius rm (see also lArmitage fc Liviol 119961 : 
iFrank et al.ll2002ll . We estima. t e n.i for X-1 using the 
expression derived bv lEggletonl (|1983ll : 

0.49ag2/3 
-(- In (1 4- ’ 

where q = Mbh/AIwr a is the binary separa¬ 
tion. Mbh is the mass of the black hole (20±4 Mq), 
and M-wr is the mass of the WR donor (26^5 Mq, 
iCrowther et al.l[2010D . Using this equation, we estimate 
rni to be 0.3-0.4a. This is con sistent with that pre¬ 
dicted by ICrowther et al.l (j20inl l. who find the Roche 
lobe radius to be ~0.4a. Although the radius of the inner 
boundary of the WR star (e.g. , where the optical d epth 
T ~ 10 ) is typically ~l -2 Rq (jLavcock et al.l[2?)T5h . the 
large outflowing winds can ex tend the “photos phere” out 
to ^8 Rq (see the review bv iCrowthen 120071 and refer¬ 
ences therein), comparable to our estimated tli. The 
outer portion of the WR “photosphere” is therefore ca¬ 
pable of fueling the BH accretion disk via Roche lobe 
overflow. 

From this estimate of the Roche lobe equivalent radius, 
we can now estimate the radius of the corona. Using a 
period P of 33 hours, an ingress time At of 8 hours, and 
an accretion disk radius of td = 0.8rLi = 0.35a, we es¬ 
timate the radius of the corona to be ~0.2a. This value 
is consist ent with corona radii of other eclipsi ng sources 
found bv IChurch fc Balucihska-Churchl (I2004D . and pro¬ 
vides additional circumstantial evidence (from the X-1 
light curve alone, and not from our spectral modeling) 
that X-1 hosts an extended corona. 

5.2. Mass Loss from Stellar Winds and Accretion onto 
the Black Hole 

We may shed additional light on the donor star by com¬ 
paring the accretion rate of the BH to the stellar mass 
loss rates. Using the unabsorbed, 0.35-8 keV luminos¬ 
ity predicted from our simultaneous eclipse egress/non¬ 
eclipsing spectral modeling ('^4.6x10^® erg s“^), we es¬ 
timate that the X-1 BH (with an inferred mass of 20±4 
Mq) is accreting at ~18% Eddington (with a lower BH 
mass radiating at a higher fraction of the Eddington 
limit; see T5.3p . A lower BH would be required to ex¬ 
plained the observed X-ray emission if X-1 possessed a 
mild degree of beaming (given the system inclination of 
z ~ 60 — 75°; iKing et ^1200111 . For example, a BH mass 
of '^10 Mq would be required to explain the observed 
X-ray luminosity with a beaming factor b ~0.5. 

The mass accretion rate M, which can be estimated 
from Lx = rjXIc^ (zy, the radiation efficiency, is assumed 
to be 0.1), is ^8x10“® Mq yr“^. Since our spectral 
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TABLE 6 

Simultaneous Fitting of ObsID 16028 Eclipse Egress and Non-Eclipsing Spectra 


Parameter 

Eclipse Egress 

Non-Eclipsing 

(1) 

(2) 

(3) 

(cm 

4.09x10^*^ (Galactic; flxed) 

fcTin 

l.libO.9 IteV 

photon temperature 

0.1 IieV 

electron temperature 

40 IceV 

optical depth 

0.17 

X^/dof 

69/70 

unabs. 0.35-8 flux 

9.6x10“^^ erg cm“^ 

unabs. 0.35-8 luminosity 

4.6x10^® erg 8“^ 

partial covering % 

86l“ 

44lj0 

-^H,cov 

12.312.6 10^^ 

1.411K ^ 10^^ 


modeling revealed no significant changes in unabsorbed 
flux, we can assume that there are no significant changes 
in the mass accretion rate, and that the observed change 
in flux is due almost entirely to the presence of high 
optical depth obscuring material during the eclipse. 

We can additionally estimate the wind mass-loss 
rate of the WR star , towr, using the approach of 
iBelczvnski et al.l (|2013ll : 

/ r X 1.5 

rhwR = 10"^^ \ T7) 

ICrowther et al.l (|2010[) measure log(L/L0) = 5.92 for 
the WR star; inputting this value into the above equa¬ 
tion yields a wind mass loss rate of 7.6 x 10“® Mq yr“^. 
This wind mass loss rate is based on mo re detailed stel¬ 
lar ev olution models than were used in ICrowther et al.l 
(|2010r) . with a comprehensive treatment of WR winds. 
As a result, our estimated mwu is roughly a fact or of 
ten higher than predicted bv ICrowther et ahl 1)201(1( 1 and 
is nearly a thousand times larger than the BH mass ac¬ 
cretion rate M. Thus, the WR winds are easily sufficient 
to fuel the BH accretion disk. 

If the donor star were not the WR star, but instead one 
of the other stars observed within the X-ray error circle of 
our HST imaging, the above equation may not be valid. 
We therefore additionally estimate the wind mass loss 
rates for an AGB star and B-type main sequence star 
(both of which were found in the vicinity of the X-ray 
source in m- 

The magnitudes of the bright, red star detected within 
the X-ray error circle in our HST imaging are consis¬ 
tent with those of an AGB star. Mass loss during 
the super-wind phase of thermally pulsating AGB stars 
(TP-AGB) can reach up to 10“® — 10“"^ Mq yr“^ (see 
iRosenfield et al.ll2ni4 and references therein), compara¬ 
ble to mwR calculated above, wi th wind speeds on th e 
order of a few tens of km s—1 ((Mattsson et al.l 1201011 . 
Thus, an AGB wind could power the observed X-ray 
emission. However, if X-1 did host an AGB donor, it 
would additionally imply a significantly lower BH mass. 

The mass loss rates for several OB mai n sequence 
stars in the Spitzer Gygnus X Legacy Survey (iHora et al.1 
[2?inl were estinrated to be a few 10 ® Mq yr ^ 
(jKobulnickv et idl 120101 1 . Although significantly lower 
than our estimated WR wind rate, only ~8% of the OB 
stellar wind would be required to fuel the X-1 BH. How¬ 
ever, the orbital period of X-1 makes this scenario un¬ 
likely. Using Equation [H we can estimate tli of the 


donor if the masses of t he binary components are known. 
The X-1 mass function (jCrowther et al.ir201(lll and an as¬ 
sumed donor mass of ~10 Mq (for typical main sequence 
B star) implies a lower black hole mass of ^6.5 Mq. 
These two masses combined with the observed orbital 
period imply a binary separation of ^ 13i?0. With this 
orbital separation, it is unlikely the main sequence star 
would fill its Roche lobe, and one would expect the BH 
to be wind-fed and not accreting via Roche lobe overflow. 
Such wind-fed systems are typically two orde r s of mag¬ 
nitud e fainter than X-1 (iLiu &: Bregmanll2005t IlIu et al.1 

[IMl). 

5.3. The Origin of the He II X4686 Emission Feature 

If X-1 is truly a BH accreting via Roche lobe overflow 
from a WR donor, the orbital separation of the binary 
is small enough that the WR star would be very close 
to filling its Roche surface. Focu sed stellar wind inodels 
of Robe lobe-filling donor stars ((Friend fc CastoH 119821 : 
iGies fc BoltonlflO^ predict that the stellar winds will 
be highly asymmetric, with the greatest mass loss occur¬ 
ring in the direction of the BH. One of the predictions of 
this model is that the He II A4686 emission feature will 
be present in the optical spectrum and display radial ve¬ 
locity variations. While it was originally thought that 
such an emission feature could only be produced in the 
densest portions of very massive blue supergiant winds 
(and not the stellar photosphere), this feature may also 
originate near the inner La grangian point where the stel¬ 
lar winds are most dense (^Walbornl]l971l : iGonti fc Feed 
11974 liHein fc CastoHll97f^ 

The He II A4686 emission feature was observed in the 
optic al spectrum of the X-1 WR donor ((Crowther et al.l 
[201(11 . along with radial velocity variations that are qual¬ 
itatively similar to those predicted by a focused wind 
through a Roche lobe over flowing donor s t ar. W e use 
the information provided in ICrowther et all ((201(111 to re¬ 
calculate the phase of each observation, assuming an or¬ 
bital period of 33.0 hours with phase 0 corresponding 
to the start of our Chandra Obs ID 16028 observation. 
Figure (TU] shows our results: high radial velocities cor¬ 
respond to the deepest part of the X-ray eclipse, with a 
minimum velocity occurring near the location of superior 
conjunction. Thus, the velocities are tracing either the 
excited wind in the vicinity of the BH or winds from the 
BH accretion disk, and not the motion of the WR star. 
This finding suggests that the BH mass measurement, 
which assumes the He II A4686 emission line is providing 
the velocity of the WR star, may not be reliable. This is 
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somewhat different from IC 10 X-1, where a phase shift 
of 7r/2 rad between the He II A4686 emission line and the 
X-ray light curve has been interpreted as originating in a 
persist ent spiral shock stru cture in the winds of the WR 
donor (iLavcock et al.l[20T5l l. 
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Fig. 10.— The He II A4686 velocity measured bv ICrowther et aTI 
ii2oTg'i . with phase corrected such that phase 0 corresponds to the 
start of our Chandra Obs 16028 observation and the assumed pe¬ 
riod is 33.0 hours. 

Our phasing of the X-ray light curve with the He II 
A4686 velocities suggest that, as the BH moves to oppo¬ 
sition, it is moving through a progressively higher col¬ 
umn density of the stellar wind. The dense portions of 
the wind primarily obscure the soft X-rays being emit¬ 
ted by the corona (which, due to its extended nature, 
has access to the cool outer portions of the accretion 
disk); a grazing eclipse by the WR donor may or may 
not also be occurring. As the BH moves from opposition 
to conjunction, the column density of obscuring material 
decreases and less of the corona and accretion disk are 
covered from the perspective of the observer. A similar 
toy model of how a focused wind from a Roche lobe over¬ 
flowing donor can create such a change in column density 
w as presented for the B H-HMXB Cyg X-1 (see Figure 8 
of iGrinberg et al.lI2015L and references therein), which 
contains an early supergiant companion (|Walbornlll973[l 
in an even tighter bina ry orbit than the case of X-1 (~5.5 
days, iGies et H I 1200811 and at a lower inclination angle 
{i ^ 27°: lOrosz et al.ll2011h . The observation that the X- 
1 He II A4686 radial velocity shows maxima and minima 
at phases where the BH is expected to have the smallest 
radial velocities further supports this picture: at opposi¬ 
tion, the dense winds are moving maximally away from 
the observer, and at conjunction the winds are moving 
maximally towards the observer. 

This model is, of course, overly simplistic. First, the 
above description does not include the effects of clumps 
in the WR wind, which may significantly contribute to 
the X-ray variability. The close proximity of the X-ray 
source will additionally alter the wind structure dramat¬ 
ically, as a large portion of the wind is expected to be 
highly ionized. The wind velocity is additionally affected 
by ionization, as the winds of WR and other massive 
stars are line-driven (iGastor e t al.lll9 75ll. which in turn 
affects the wind density (|Hatchett fc McCravlll977ll . An 
asymmetric wind was inferred for Cyg X-3 from radial 
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velo city measurements o f the He I 2.0587^m absorption 
line (|Hanson et al.l[2000ll : however, the He I line obtained 
in the X-1 optical spectrum was too we ak for a similarly 
detailed analysis (jCrowther et al.l 1201(1(1 . Detailed mod¬ 
eling of the WR wind structure is beyond the scope of 
this work, and the limited X-ray observations available 
for X-1 make it unlikely that we would be successful in 
distinguishing bet ween different models. 

It was noted by iCrowther et al.l (|201Clll that the He 11 
A4686 equivalent width (W\ ~ 56A) is a factor of ~2 
lower than for similar WR stars in the Milky Way. This 
low value of the equivalent width is possibly due to dilu¬ 
tion from other (unresolved) stars along the line of sight. 
The cool giant star within the X-ray error circle may be 
the primary contributor to this dilution, as such stars are 
not hot enough to produce He II A4686 spectral features. 

6 . SUMMARY 

We have presented new Chandra and HST observations 
of the WR -I- BH binary NGC 300 X-1. The periodic 
dips in the X-ray light curves has led previous authors 
to conclude that the X-ray source experiences a grazing 
eclipse by the WR donor. We capture a large portion of 
an eclipse egress in one of our new Chandra observations, 
and find that intervals of low X-ray flux correspond to a 
harder underlying spectrum. Simultaneous spectral fit¬ 
ting of the eclipse egress spectrum and a portion of the 
spectrum corresponding to superior conjunction of the 
BH are consistent with a partially covered multicolor ac¬ 
cretion disk and an extended, Comptonized corona. Dur¬ 
ing the eclipse egress, X-1 is partially obscured, possibly 
by a cloud in the donor wind, with a neutral hydrogen 
equivalent column density of ~10^^ cm“^ and a partial 
covering fraction of ^86%; however, this absorbing col¬ 
umn density decreases by roughly an order of magnitude 
and the covering fraction is reduced by half by the time 
the X-ray source has reached superior conjunction. We 
argue that the observed change in the column density 
and covering fraction with orbital phase is due to the BH 
moving through dense clouds in the wind of the donor 
star. 

Our HST observations reveal many optical sources 
within the X-ray error circle, including the WR star, an 
AGB star, and what is likely a high mass main sequence 
star. We consider the possibility that the donor star in 
the X-1 system is an AGB star or high mass main se¬ 
quence star. While we can rule out the possibility of a 
main sequence donor, we cannot rule out the possibility 
that the AGB star is the true donor to the BH. Our obser¬ 
vations suggest that NGC 300 X-1 is either a high-mass 
WR -I- BH binary or a low-mass AGB -|- BH binary, likely 
with an extended corona, where the observed changes in 
X-ray flux are largely due to the BH’s motion through 
the dense, clumpy stellar wind of the donor star. 

The He II A4686 emission feature and radial velocity 
varia tions observed in the optical spectrum of the WR 
star (jCrowther et al.l[2010ll are similar to those predicted 
by focused wind models. When we rephase the He II 
A4686 radial velocity measurements to match our folded 
X-ray light curve, we find the highest radial velocities 
correspond to the deepest part of the X-ray eclipse, with 
minimum velocities occurring near the location of supe¬ 
rior conjunction. This is consistent with the emission fea¬ 
ture originated from the dense material flowing through 
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the inner Lagrangian point or originating from the BH 
accretion disk, and not in the low-density photosphere of 
the WR star. This finding suggests that earlier BH mass 
measurements, which were inferred from the He II A4686 
emission, may not be reliable. 
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